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Tb ■ Abstract 

O . Motivated by interest in pedestrian traffic we study two lanes (one- 

dimensional lattices) of length L that intersect at a single site. Each 
lane is modeled by a TASEP (Totally Asymmetric Exclusion Process) . 

^ . The particles enter and leave lane a (where a = 1,2) with probabil- 

^ [ ities aa and /3a, respectively. We employ the 'frozen shuffie' update 

introduced in earlier work [C. Appert-Rolland et al, J. Stat. Mech. 



ly-^ ■ (2011) P07009], in which the particle positions are updated in a fixed 



random order. We find analytically that each lane may be in a 'free 
C3 I flow' or in a 'jammed' state. Hence the phase diagram in the domain 

< ai, a2 < 1 consists of four regions with boundaries depending on 

Pi and (32- The regions meet in a single point on the diagonal of the 

domain. Our analytical predictions for the phase boundaries as well 

^ ' as for the currents and densities in each phase are confirmed by Monte 

Carlo simulations. 
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1 Introduction 

Pedestrian motion has raised increasing interest in the past years, 
both from a practical and a theoretical point of view. Understand- 
ing the behavior of crowds or of waiting lines is still a challenge. 
Simplified models may help to understand the behavior of individuals and 
the resulting collective behavior in various settings. In a large class of models 
[H |2] pedestrians are represented as hard core particles moving on a lattice 
according to certain rules of motion. One important ingredient of these rules 
is the type of update scheme that is employed. Actually the update scheme 
is an integral part of the model definition; changing the scheme may change 
the interpretation and the properties of the model [3j . 

In the past two types of updates have been used for pedestrians modeling: 
the random shuffle update [U |5l |6l [7] which has been later replaced by the 
parallel update [HI |9]. In [101 [H] we have proposed a new update scheme for 
pedestrian modeling, the frozen shuffle update, that we shall use in this paper. 
Its characteristic feature is that during each time step all particles present in 
the system are updated in a fixed random sequence. Newly entering particles 
are inserted in this updating sequence and exiting particles are deleted from 
it according to a suitable algorithm. Frozen shuffle update was inspired 
originally by the need for a physically motivated rule of priority in cases 
where more than one particle attempts to hop simultaneously towards the 
same target site. This update has the additional advantages that it is easily 
implemented in a Monte Carlo simulation and lends itself well to analytic 
study. 

The consequences of frozen shuffle update were worked out previously for 
the case of a one- dimensional totally asymmetric exclusion process (TASEP) 
both on a ring [TU] and with open boundary conditions [TT] . On a finite one- 
dimensional lattice with open boundaries two parameters a and /? describe 
the probabilities for particles to enter the system at one end and to leave it 
at the other end. In this case the particle density p and the current J must 
be determined as a function of a and /3. For varying a there appears to be 
a critical point a = /3 between a 'free fiow' and a 'jammed' state. 

One of our purposes is to model pedestrian motion at the intersection of 
two corridors or two streets and to study how global structures emerge from 
local interactions. As a step toward this goal we study in the present paper 
a TASEP on two perpendicular traffic lanes, 1 and 2, that intersect at a 
single lattice site and whose entrance and exit parameters are ai,a2, (3i, and 
/32- The main question is again to determine the stationary state currents Ji 
and J2 in this two lane system as a function of these four parameters. For 
each lane one may expect two possibilities, a free fiow (F) or a jammed (J) 
state. We will study the phase diagram in the aia2 plane for < ai, 0:2 < 1, 



considering /3i and P2 as fixed parameters. Indeed we find a division of this 
square domain into four different regions denoted FF, FJ, JF, and JJ, and 
separated by pliase boundaries for wliich we obtain analytic expressions. 

In earlier analytic work 10] on the frozen shuffle update the concept of 
a 'platoon' was introduceqll. It will again play a role in this paper. We 
will, moreover, point out here a new phenomenon, to be called the 'pairing 
mechanism', which is operative at the intersection. It says, basically, that 
when both lanes are in the jammed state, a platoon crossing the intersection 
on lane 1 is always accompanied by a platoon crossing the intersection on 
lane 2. This mechanism, which is an unintended consequence of the rules of 
motion, will enable us to extend the theoretical analysis from the single lane 
to the case of two intersecting lanes. 

This article is organized as follows. In section [2] we define the exact rules 
of motion of this intersecting lane model and recall the concept of a 'platoon'. 
In a short section [3] we recall the single lane results that will be needed again 
here. In subsection 14. II we argue that the FF phase, expected to exist at low 
entrance rates, cannot extend beyond a certain curve in the aia2 plane. In 
subsection 14.21 we show that if a J J phase exists, the exiting flow must obey a 
pairing mechanism. Exploiting this mechanism we determine in subsections 
14.31 and 14.41 the phase boundaries of the intermediate F J/ JF states with the 
JJ and FF states, respectively, and thereby confirm the existence of all four 
possible phases. We obtain analytical expressions for the currents in both 
lanes in each of the four phases. In subsection 14. 51 we consider various limits in 
the aia2 domain. In section [5] we derive expressions for the particle densities 
which, in contrast to the currents, are discontinuous at the phase boundaries. 
In section [6] we present a few simulation results. The data for the current 
fall right onto the theoretical curves, whereas the density data show finite 
size effects similar to those encountered and explained in the single lane case 
[llj. Section [7] is our conclusion. 

2 Rules of the motion 

2.1 Rules 

We consider the geometry shown in figure [H consisting of two perpendicular 
one-dimensional lattices (or: lanes) labeled by an index a = 1,2. Hard core 
particles may move to the right on lane 1 and upward on lane 2. The lanes 
have Li and L2 sites, respectively, plus a common intersection site. Particles 
enter at the two extremal sites and exit when leaving the intersection site. 
In our analytical treatment we shall imagine that both Lo- tend to infinity. 



^This term has been borrowed from road traffic. 
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Figure 1: The two lane configuration studied in this work. The arrows indicate 
the flow direction. A heavy (red) bar indicates an end-of-platoon. The entrance 
probabilities are ai and 02, the exit probabilities from the intersection site are /3i 
and f32- Inset: the single lane with parameters a and /3 studied in reference [llj . 



As usual in critical phenomena, finite size effects will be important only in 
the critical region. In practice, our infinite system results apply as soon as 
the Lfj are large compared to the boundary layers near the entrance or exit. 

A particle L when entering the system, is assigned, in a way discussed 
below, a phasci ^i € [0, 1) which it keeps as a fixed attribute until it exits 
the system. Although the time t is continuous, the time evolution is best 
described in terms of integer time steps s = 1,2,3,.... During each time 
step the particles are visited in the order of increasing phaseqj and their 
positions are updated according to the following rules. 

General rule. When a lane 1 (lane 2) particle is updated, it moves one 
lattice distance to the right (upward) if the target site is emptja, and does 
not move if the target site is occupied. The update of particle i during the 
sth time step is considered to take place at the exact moment t = s + Tj . 

The general rule must be supplemented by two special rules for entering 
and exiting particles. 

Exiting rule. A particle i which at the beginning of the sth time step is 
on the intersection site, will at time s + Xj leave the system with probability 
Pi (or P2) according to whether it has arrived through lane 1 (or lane 2). 
Once it leaves the intersection site, we do not consider it any longer: it has 
left the systeno. 

Entering rule. When the entrance site of lane a becomes vacant at time 
t, it will be occupied by a new particle, injected from outside, at a random 
time t' = t + T. Here T is drawn from the exponential distribution 

P^{T)=a„Q-'^-T^ 0<T<oo, (2.1) 

where the 'entrance rate' Oo- > is a model parameter. 

Equivalent to a^ is the 'entrance probability' a„ defined by 

a^ = l-e~"% c^ = l,2, (2.2) 

which is the conditional probability that the entrance site of lane a is occupied 
at time t + 1 given that it was vacant at time t. Henceforth we will sometimes 
let Oo- and a^ occur in the same expression. 



^We use in this work the term 'phase' both to designate the r^ assigned to the particles 
and to refer to the different types of stationary states of the system as a whole. No 
confusion need arise. 

•^For a closed system with a fixed number N of particles a randomly chosen permutation 
of the particles may replace the assignment of phases \Q. 

^The TASEP considered in this paper is deterministic in the bulk. For particles hopping 
forward to an empty target site with a probability p < 1, analytical predictions would 
probably be more difficult. 

^In fact, it has only left our window of observation: we could consider that each lane ex- 
tends beyond the point of intersection, and that the particle, once it leaves the intersection, 
enters a free flow state where it continues to move at every time step. 



To see the motivation for the above entering rule, one may notice that each 
particle arrival on the entrance site (or, for that matter, on any other site) is 
followed by one unit of 'dead time' during which no new arrival on that site 
is possible. Subject to this dead time condition, the entering rule distributes 
the instants of arrival of the particles at the entrance site uniformly on the 
time axislj. 

2.2 Platoons 

The phases Tj may be regarded as quenched random variables. With the rules 
stated above the particle motion is deterministic apart from the stochastic 
phase assignment at the entrances and, when /3i < 1 or /32 < 1, the stochastic 
exits. As a consequence of the entering rule the phase r' of the newly injected 
particle is related to the phase r of its predecessor in the same lane by 

r' = (r + T) modi. (2.3) 

It follows that there are correlational between the phases of successive par- 
ticles in the same lane. Following a lane in the direction opposite to the 
particle flow, one may group the particles together into sequences of increas- 
ing phases, a phase decrease signaling the beginning of a new sequence. When 
particles constituting such an increasing phase sequence occupy consecutive 
sites, they are said to constitute a 'platoon'. The average platoon length u 
associated with an entrance probability a is given by [11] 

Til . ^ 

1 + . (2.4) 



//(a) a a 

where a = 1 — e"''. This quantity will play an essential role in the analysis 
that follows. 



3 Stationary states in a single lane 

The single lane problem with boundary conditions a and /3, shown in the 
inset of figure [H has yielded [H] results some of which will again be needed 
here. First, we know that the entrance probability a (for large enough /3) 
imposes a 'free flow' bulk statqj, that is, one in which all attempted moves 
are successful, which carries a current 

J^(«) = -^. (3.1) 

1 + a 



^An analogy with a one-dimensional system of hard rods of unit length was pointed 
out in reference |11| . 

^Described in detail in reference [TT]. 

® Except for a jammed boundary layer of fluctuating size near the exit. 
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Secondly, it was shown yjLj that the exit probabihty (3 (for large enough a) 
imposes a jammed bulk statqj. This is a state in which all particles belong to 
platoons and successive platoons are separated by at most a single vacancy. 
The jammed state has an outgoing current 

jJ(a,/3) = ^, (3.2) 

with z/ determined by a through (12. 4p . This exit flow can be sustained if 
and only if the entering flow is sufficiently large, that is for J^ > J\ The 
equality J^(a) = J'^(a,/3) therefore defines a critical point, which turns out 
to occur for a = p. As a consequence, the stationary state current J is equal 
to J = J^{a) ioT a < (3 and J = J'^{a,/3) for a > /3; at the critical point 
it is continuous but undergoes a change of slope. Finally, for a = /3 the two 
phases coexist in the system and are spatially separated by a sharp domain 
wall. 



4 Phase diagram of the two lane system 

4.1 The FF phase 

In the two lane system the entrance probabilities ai and 0:2 strive to impose 
independent free fiow states in each lane, that is, an FF phase with currentdlj 



1 + a^ 



a = 1,2. (4.1) 



The two currents interact at the intersection site where moreover they are 
subject to random exits with probabilities /3i and /32. We anticipate that if at 
given /3i and /32 the entrance probabilities ai and 02 become small enough, 
the system will be in an FF phase. The interaction between the currents at 
the exit site may then occasionally delay individual particles, but will not 
create waiting queues that grow without limit. 

However, the rules of the motion are such that at each time step at most 
a single particle can leave the exit site. This immediately yields a bound for 
the FF phase in the aia2 plane: whenever J^(ai) + J^{a2) > 1, there must 
necessarily occur formation of an ever growing waiting line in at least one of 



^Except for a free flow boundary layer of fluctuating size near the entrance. Only in 
the limit L — >■ 00 are the two phases well-defined in the sense that tunneling between them 
becomes impossible. 

^'^A symbol with a single upper index, F or J, refers to an auxiliary one lane system; 
a symbol with a double upper index refers to one of the lanes a = 1 , 2 of the two lane 
system under study. 
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Figure 2: Sequence of snapshots illustrating the pairing effect of platoons at the 
intersection site. In this example ni = 3 and n2 = mi = m,2 = 2. See text. 

the two lanes and the system cannot then be in its FF phase. This condition 
may be rewritten as aia2 > 1. In subsection 14.41 we will show by explicit 
calculation that the FF phase does indeed exist and analytically determine 
its phase boundary. 

4.2 Pairing mechanism 

There is no standard way of finding the phase diagram of this two lane 
system. We therefore develop following reasoning. 

Let us suppose now that both lanes are jammed, that is, the system is in 
a JJ state. In the two-lane problem there then appears a new phenomenon. 
The rules of the motion have an unintended consequence that we will call 
the pairing mechanism. This is the phenomenon that the exiting platoons 
of the two lanes are rigorously paired: for each platoon exiting from lane 1 
there is also a platoon exiting from lane 2, and vice versa. To demonstrate 
this effect we refer to figure O 

Figure [2^ shows a particle configuration near the exit at some integer 
instant of time t = s at which the intersection site is empty. A heavy (red) 
bar (below or to the left of a particle) marks an end-of-platoon. In each 
lane there is a platoon (dark colored particles) waiting to enter the empty 
intersection site; the platoons have lengths Ui and ^2 and are headed by the 
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particles marked 1 and 1'. At the next time step, t = s + 1, the platoon head 
with the lower phase will hop onto the intersection site, its whole platoon will 
advance by one lattice distance, and it will block the other waiting particle. 
Let us suppose, as shown in figure [2)d, that it is the horizontally moving 
platoon that advances. During each of the next time steps, the lane 1 particle 
(marked 1), as long as it occupies the exit site, will exit with probability /3i, 
so it will exit only after on average /3f^ time steps. During the same time 
step in which it exits, out of the two particles (marked 2 and 1') that are 
waiting to hop onto the exit site, again the one with the lower phase will 
effectively hop, pull along its whole platoon, and block the other one. Let 
us suppose that it is the lane 2 particle (marked 1') that advances. This 
takes us to the configuration of figure [2t. It is equivalent to the one of [2)3, 
except for an interchange of the roles of lanes 1 and 2, accompanied by the 
replacements rii i— )■ n2 and n2 H- ni — 1. The procedure taking us form[2}3 to 
^ will now repeat itself mutatis mutandis and each time either ni or n2 will 
decrease by one unit. Lane 1 and lane 2 particles have average exit times 
/3]"^ and /3^\ respectively. At some point the last particle of one of the two 
platoons is on the intersection site. Let us suppose this is a lane 2 particle, 
as in [2ji where it is marked 2'. During the same the time step in which 2' 
leaves the intersection site, that site will be occupied by the particle waiting 
in the other lane (marked 2), which has a higher phase. This will result in 
the situation of figure [2^. The next lane 2 particle (marked 3') belongs to 
the next platoon; if it has already arrived at the waiting site (which may or 
may not be the case), it will be blocked in that time step. It will similarly be 
blocked in all following time steps, until the last remaining particle (marked 
3) of the lane 1 platoon leaves the intersection site. The situation that then 
results is depicted in figure [2]F. It is identical to that of figure [2^, except that 
now the next two platoons, of lengths mi and 777-2 and headed by particles 4 
and 3', are waiting to enter the intersection site. This is the pairing effect. 

We remark parenthetically that this pairing argument is easily extended 
to an arbitrary number p of lanes intersecting at a single site, when they are 
all in the jammed phase. We will not, however, attempt to consider here 
such more general geometries. 

We will now exploit this effect to find an expression for the current in the 
J J phase. In order to arrive at the situation of figure [2]F starting from the one 
of figure [2^ there is first the time step in which the exit site gets occupied. 
Next, there are 77i lane 1 particles and 772 lane 2 particles that leave the 
intersection site subject to the exit probabilities /3i and /32, respectively. The 
total time t„j„2 needed for this process and averaged over all exit histories 
therefore is tnin2 = ^i/^f^ + '^-2/32^"'^ + 1- Let ui = z/(q;i) and z/2 = 2/(0:2) be 
the average platoon lengths in the two lanes. Then the mean exit time texit 
of an arbitrary pair of platoons to exit is the average of tnin2 over all platoon 



lengths. This yields 

texit = ^ + 7^ + 1- (4.2) 

Pi P2 

Hence in the JJ phase the outgoing currents J('^ and Jg'^ of the two lanes are 
given by 

4^^= ^/; , 1 ^ - = 1,2. (4.3) 

This expression is a nontrivial generalization of the single lane formula fl3.2p . 
Both currents fl4.3|) depend on all four parameters ai, 0^2, /3i, /32. The ratios 
Ufj/P^, which will reappear frequently below, show the scaling with P„ of the 
time that the intersection site is occupied by lane a particles. 

All elements are in place now for us to go on and find the phase boundaries 
in the domain < ai,a2 < I- 

4.3 Boundaries of the JJ and FJ/JF phases 

Let us suppose the system is in the JJ phase. The condition for the system 
to be able to sustain these queues is that in both lanes the out-current J^'^ 
be smaller than the corresponding free flow entrance driven current J^{aa)- 
That is, for the JJ phase to be stable we should satisfy the two inequalities 

J^(«.)>Ji^ o- = l,2, (4.4) 

or, upon substituting (14.11) and (14. 3 p in (14.40 . 

'^ > . , t , . , - = 1,2. (4.5) 



On the borderline of the JJ phase equation (14. 5 p should hold as an equality 
for either a = 1 or a = 2. To rewrite this equality we invert both of its 
members and use (12. 4p . It then follows that 

+ ^, f^=l,2, (4.6) 



Oia Pi P2 



or, equivalently, 



"■^(-4) = w- (^■''') 

V"2 P2J Pi 

Equations (14. 7p represent two intersecting curves in the a 102 plane. Although 
they depend on the parameters /3i and /32, their point of intersection always 
lies on the diagonal ai = 02. To see this, note that z/g. = vi^aa) is a function 
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Figure 3: Phase diagram in the 0102 plane for /3i = /32 = 1- The heavy sohd 
Hnes are phase boundaries. Dashed hne: the diagonal. The symbols F (free) and 
J (jammed) refer to lanes 1 and 2 in the order given. The four-phase point is at 

(01,02) = (|,|). 



only of Ua and hence a\ = a2 = Oc implies that v\ = V2 = Vc- Using this in 
(14.71) we see that v^. divides out and that both equations are satisfied by 



(Xr 



/3i/32 
/3i + /32 



(4. 



Hence equation (I4.7ap gives the JJ/FJ boundary in the triangle above the 
diagonal a\ = 02 and fl4.7bp gives the JJ/JF boundary in the triangle below 
this diagonal. For the special case f3i = (32 = I these phase boundaries 
are shown in figure Ej which is symmetric with respect to the diagonal. An 
example of the general case with /3i 7^ P2 is shown in figure HI where this 
symmetry has been lost. 



4.4 Boundaries of the FF and FJ/JF phases 

We wish to find now the borderline between the FJ and JF phases, on the 
one hand, and the FF phase on the other hand. To be definite, let us suppose 
the system is in the FJ phase so that we know that the current in lane 1 is 
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Figure 4: Phase diagram in the a\a2 plane for /3i = 0.95 and /32 = 0.45. All 
symbols are as in figure O The four-phase point is on the diagonal at {pcc-, ol^ with 
Uc given by (j4.8p . As shown in section 14.5.2^ the boundary curves of the FF phase 
intersect the a\ and 02 axes in a\ = (3i and 02 = (32, respectively. Dotted circle: 
path along which the simulation results of figures [6] and [7] were obtained. 
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given by its free flow expression 

Jf J = J^(ai). (4.9) 

Tlie current Jj^^ in lane 2 is, however, unknown. In order to calculate Jj'^ 
we cannot invoke now the pairing mechanism, since it requires both lanes 
to be jammed. Instead, let us suppose that for every platoon that exits the 
jammed lane 2, and that is known to contain on average z/2 walkers, there are 
on average /ii walkers that exit lane 1. Here fii is unknown but necessarily 
satisfiecj 

/ii < ^^1 . (4.10) 

Then we have for the exiting currents in the FJ phase the two asymmetric 
expressions 

/3i "^ fe "^ "- 
•'" = ITT^- (''■lib) 

/3l "^ /32 "^ ^ 

Upon combining (14. 9^ with (14.1 lap and using (13. ip we may solve for /xi and 
find 



The condition for the sustainability of the jammed phase in lane 2 is 

J^'Myj!'. (4.13) 

The phase transition line is obtained when (I4.13P holds as an equality, which. 



with the substitutions of (14.1 Ibp and ( 13. ip . happens for 



02 1^2 



l+«2 t+t + 1 



1^2 



a, 1+ i-ir)«i 



1 + 0-1 



(f + «1 



1^2(^2 '^['-T^ 



1+ l-±]a. 



(4.14) 



1^2 + (^2 1 + Oi 

""^""^ Relation (|4.10p . valid in the FJ phase, becomes an equality when lane 1 also gets 
jammed, so that fj,i — vi should give again the JJ/FJ boundary. This may be verified by 
explicit calculation. 
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where to pass from the first to the second hne we first noticed that by (14. 11 aft 
and (14. 9 p the denominator on the RHS is equal to /ii/Jf"' = ni/ J^{ai) = 
/ii(l + ai)/ai and then substituted for /ii expression (I4.12p . We may solve 
(I4.14P for «! and find 

"^-l-/3i(l-i?.)' ^'-''^ 

where we introduced the abbreviation i?2 = -R(a2,/^2) with 

We will continue to consider /3i and P2 as fixed parameters. When in (14.150 
we use (I4.16P for R2, (12. 4p for 1^2, and (12.21) for 0:2, it becomes an exphcit 
solution for ai in terms of 02, or equivalently, for ai in terms of 0:2 . Hence 
(I4.15P constitutes our final result for the FF/FJ boundary. A permutation of 
indices gives the FF/JF boundary. These phase boundaries are again shown 
in figures [3] and H] for the symmetric case with /3i = (32 = I and for a typical 
asymmetric case, respectively. 

4.5 Limiting cases 

We consider in this section the limit behavior of the phase boundaries as 
they approach the borders of the domain < ai, 0^2 < 1. 

4.5.1 Boundaries of the JJ phase 

One obtains from (I4.7al) the behavior of the JJ/FJ boundary in the limit of 
small «! by noticing that limQ,j_!.o i^i = 2 and that therefore in that limit the 
LHS of (I4.7ap diverges, which forces z/2 on the RHS also to diverge. Using 
next that for 02 — >■ 1 one has 1^2 — — log(l — 02) one finds 

02 ^ 1 - e-2^2/"S ai -^ 0. (4.17) 

A permutation of indices gives the asymptotic behavior of the JJ/ JF bound- 
ary in the limit 02 ~^ 0. The exponentials of the inverse functions 1/ai and 
1/0:2 explain the extremely rapid alignment of these curves along the edges 
of the figure. 

4.5.2 Boundaries of the FF phase 

We wish to find the point of intersection of the FF/JF (FF/FJ) boundary 
with the horizontal (vertical) axis. It is located at ai = /3i (at 02 = ^2)- To 
show this, we ask what the limiting value of q;2 is when ai -^ 0. It is useful 
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to notice that the quantity R2 that occurs in fl4.15p may be expressed as a 
ratio of two single lane currents, 

R. ^ f:^^i , ■ (4.18) 

The 'physical' argument goes as follows. For ai = lane 1 is unoccupied 
and the intersecting lane problem reduces to that of the single open-ended 
lane with boundary conditions (02? /32); whose critical point is known [11] 
to occur at 02 = ^2- Mathematically, ai = implies ai = 0; when this is 
substituted for the LHS of f l4.15p we find R2 = 1, after which f l4.16p yields 
J^^'^'^{a2) = J^^^{(y2, (32)- When this equality is worked out we obtain the 
same result ^2 = /32- Finding the limit behavior for a2 ^ amounts to a 
permutation of indices. 

The straight line (not drawn in figures [3] and H]) that connects these two 
points of intersection has the equation 

and also passes through the critical point (ac,«c)- In both figures, |3]and|ll 
the boundary delimiting the FF phase is shghtly curved and falls just below 
this straight line. 

5 Particle density 

The determination of the phase diagram was based exclusively on the analysis 
of the particle currents in the different phases. From the preceding construc- 
tion it follows that the currents are continuous at the phase transition lines. 
This differentiates them from the particle densities, which are the quantities 
of interest in this section. We will denote a particle density generically by 
the symbol p, to which we attach indices according to the same convention as 
used for J. In all phases we have the relation J = vp, where p is the particle 
density and v the average particle velocity. Since in the free flow phase all 
particles have velocity f = 1, we have p^ = J^ and therefore 

(5.1) 



pf 


"cr 


a = 


1,2, 


1 + a, ' 


I^J 


ai 


Pf = - 


^2 


'1 — 


1 , „ ' 


1 1 „ 



(5.2) 
1 -|- ai i -t- U2 

In the jammed phase we have generically the relatioio J = ^(1 — p), where 
p is as before the average platoon length. This gives the three relations 

4' = u^{l-pf), a = 1,2, (5.3) 



""^^This relation, derived in [10] . is a direct consequence of the structure of the jammed 
state described in section [Sj 
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Jl' = v,{l-pl'), Jf = z.i(l-pf). (5.4) 

Solving these for the densities using the expressions found in section H] for 
the currents we get 

£L -L £2. 

P'J= .J\Jl, ^ ^ = 1,2, (5.5) 

/3i "^ /32 "^ ^ 

An _|_ £2. VI I A'2 

FJ _ h ^fe JF _ h ^fe /, ^^ 

/3i "^ /32 "^ -^ /3i "^ ft "^ -^ 

Remarkably, equation (15.51) shows that in the JJ phase the particle densities 
in the two lanes are equal irrespective of the values of cti, a2, /3i, /32. 

Equations (15.11) . (15. 2p . (15. 5p . and (l5.6p constitute our analytic results for 
the particle densities in the four different phases. 



6 Simulations 

In order to test the theory of the preceding sections we determined the phase 
of the system for fixed /3i = /32 = 0.6 on a grid of points in the aia2 plane. 
The grid was refined in a square region around the four-phase point, which 
according to equation (14. 8 p occurs at {ai,a2) = (0.3,0.3). To determine the 
phase for a specific pair (ai, 02), simulations were performed on intersecting 
lattices of lengths Li = L2 = L = 75. In a finite system the entrance and 
exit boundary conditions try to impose different phases, which as a conse- 
quence will be separated by a domain wall pT| [T2| IT3| [H] . Away from the 
critical point the fluctuating domain wall position will be localized within 
some finite distance from one of the lane ends; upon approach of critical- 
ity this localization length increases until it attains the lane length L. In 
our simulations the domain wall position was determined in each lane as de- 
scribed in [11] . We then averaged it over 5 000 000 time steps after having 
first discarded a transient of 5 000 time steps in order to make sure that the 
system was stationary. The lane was classified F or J if its mean position 
was closer to the exit or closer to the entrance, respectively. The results are 
represented in figure [5l They are in perfect agreement with the theoretical 
phase boundaries, within the resolution of the grid. 

A more detailed simulation was carried out for the asymmetric case of 
figure m In the aia2 domain we considered the circular path of radius 0.15 
and centered in (ac, oJc), represented in figure HI In each one of 64 equidistant 
points along this circle we determined the stationary state densities pi and 
P2 as well as the currents Ji and J2 in the two lanes, for lane length L = 600. 
Each data point was obtained by first discarding a transient period of 10 000 
time steps and then averaging over 100 000 time steps. The results, together 
with the theoretical predictions, are shown in figures E] and [TJ where is the 
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Figure 5: Phase diagram in the aia2 plane for /3i = /32 = 0.6. The nature of the 
phase was determined (see text) on a grid of points with a denser covering near 
the four-phase point (01,02) = (0.3,0.3). The phase boundaries are those given 
by theory. 
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Figure 6: The current, denoted by the generic symbol J, as a function of the angle 
(h in lanes 1 and 2. 



17 



n 


lane 1 


simulation 


X 


lane 2 


simulation 




lanes 1 & 2, theory | 




lane 1 


tlieory 




lane 2 


tlieory 




Figure 7: The density, denoted by the generic symbol p, as a function of the angle 
(f) in lanes 1 and 2. The black part of the theoretical curves is common to both 
lanes. Comparison to figure [6] shows that for a lane in the free flow phase we have 
J = p. 



angle between the radius vector on the circle and the positive ai axis. Dotted 
vertical lines indicate the positions of the phase boundaries, labeled by the 
same lettering A, B, C, D as in figure HI The error bars of the data points are 
smaller than the symbols. The current data of figure E] fall perfectly on the 
theoretical curve throughout the whole range of measurements. The density 
data of figure [7] show a clear deviation from the theoretical prediction at those 
phase boundaries where the latter is discontinuous. Indeed, the prediction is 
for lane length L = oo and these deviations appear to be finite size effects. We 
have verified that indeed this effect decreases when L goes up. Their physical 
origin is the same as was found for a single lane [H]: they result from the 
formation of a fiuctuating jammed boundary layer near the exit (when the 
lane is in its free fiow phase), or of a fiuctuating free fiow boundary layer 
near the entrance (when the system is in its jammed phase). 

7 Conclusion 

We have studied pedestrian traffic on two semi-infinite one-dimensional lat- 
tices, or lanes, that intersect in a common end point. The pedestrians are 
modeled as the particles of a TASEP, that is, as hard core particles capable 
of moving only in a single direction, in the present case toward the exit. 
When leaving the intersection site, a particle exits the system. The particle 
positions were updated with 'frozen shuffle' dynamics, described in section 
12.11 and argued to be a natural choice for pedestrian motion. The updat- 
ing is easy to implement in a Monte Carlo simulation an also lends itself 
particularly well to analytical study. 

Each of the lanes (labeled by a = 1, 2) is characterized by a parameter a^- 
governing the entrance of the particles at — oo, and another one, /3a, governing 
their exit from the intersection site. For arbitrary fixed f3i and /32 we have 
determined analytically the phase boundaries in the aia2 plane. It appears 
that each lane may be in either a free flow (F) or a jammed (J) phase, 
which results in a partition of the phase diagram in the aia2 plane into 
four regions, JJ, FJ, JF, and FF. Explicit expressions have been found for 
the phase boundaries between these regions. An essential element in our 
analysis is the pairing effect that we have shown to occur when both lanes 
are in the jammed phase: in that case each platoon exiting from one lane 
is accompanied by a platoon exiting from the other lane. Once the pairing 
effect is established, the analytical expressions for the macroscopic quantities 
of greatest interest, the currents and the particle densities, become accessible 
via reasonably simple mathematics. We have determined them analytically 
for each region of the phase diagram. All our analytical findings have been 
corroborated by Monte Carlo simulations presented in section O 

This work is to be seen as a first step towards the study of the intersection 
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of larger corridors, in which case the possibihty of lateral hops may also have 
to be included. We will leave the study of such more complicated geometries 
and hopping rules to future work. 
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